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For a long time a challenge for mathematical physics has been and is still the idea, due to 
Hilbert at the dawn of the quantum age, to relate the non-trivial zeros of the Riemann zeta 
function with a spectrum of a Hermitian operator in a Hilbert space. 

The Riemann zeta function ([s) is defined for complex s = a + ip and Res > 1 by the 
equation 

1 

The Riemann hypothesis [1] states that the non-trivial complex zeros of ({s) lie on the critical 
line Res = 1/2. 

This is a very good motivation for understanding the original Riemann conjecture to find a 
self-adjoint operator d (or a positive one d{l—d)) defined in a Hilbert space and with a spectrum 
given by the non-trivial zeros of the zeta function. In favour of this idea are the rigorous results 
[2] , [3] on the spectral theory of the Laplace - Beltrami operator on the Poincare complex upper 
half-plane, showing that the eigenf unctions are automorphic with respect to a discrete group 
of fraction- linear transformations. The most important application of representation theory to 
automorphic forms is the work of Selberg [4J who has shown that, if the fundamental domain 
of the discrete subgroup T C SL{2, R) is compact, then the spectrum of the elliptic operator 
for which the eigenvalue equation 

has F-invariant solutions is determined by the real values pn given by Cr(l/2 + ^p/2) = 0, where 
Cr{s) is a meromorphic function of s on the whole plane. The spectrum is continuous if the 
fundamental domain has a finite volume. 

There has been theoretical [S], [T3] and experimental [B], [7] evidence in support of the 
Riemann hypothesis, based on the fluctuations of spacing between consecutive zeros of zeta. 

In a recent paper, Connes [8j has given a spectral interpretation of the critical zeros and 
a geometrical interpretation of the explicit formulas of number theory as a trace formula in a 
non-commutative space. 

This work was inspired by the of idea of Berry and Keating, discussed in Cambridge [9j and 
in a subsequent paper [10], that the real solutions En of 

C{l/2 + tEn) = (3) 

are energy levels, eigenvalues of a quantum Hermitian operator (the Riemann operator) asso- 
ciated with the one-dimensional classical hyperbolic Hamiltonian 

Hci{x,p) = xp, (4) 

where x and p are the conjugate coordinate and momentum. They suggest a quantization 
condition generating Riemann zeros, for which however they 'see no way to interpret geomet- 
rically'. 
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We propose that such a condition can be consistently obtained as a boundary quantization 
condition of a hyperbohc dynamical system within conformal geometry. 

The characteristic properties of the classical dynamics (the Riemann dynamics) and the 
quantum analogue are listed and commented in pJOj. It has the main features of a chaotic 
dynamics, namely instability, complex periodic orbits, no time-reversal symmetry, this in turn 
makes the quantum-classical correspondence not yet completely understood. 

A prototype for the development of a theory of hyperbolic systems has been the free motion 
on a surface of constant negative curvature (a pseudosphere) , which is one of the first models 
of chaotic motion. On this surface there exists a well-defined quantum dynamics where the 
Laplace - Beltrami operator acts as the Hamiltonian operator. The pseudosphere is modelled 
by the Poincare unit disc or by the complex upper half-plane, which is the conformal image of 
the unit disc under a particular fraction-linear transformation [11]. The metric is conformal, 
proportional to the plane Euclidean metric at each point. This suggests an analysis of the 
hyperbolic dynamical system and its quantum analogue by applying the methods of conformal 
representation theory. A motivation for this is the geometry of points at infinity (or equivalently 
all directions along the cone), namely the property that any isometry of the hyperbolic space 
can be extended to a conformal diffeormorphism of the boundary sphere, which is known as 
the 'Moebius transform'. 

The classical Hamiltonian H = xp of linear dilation, i.e. multiplication in x and contraction 
in p, gives the Hamiltonian equations: 

X = X, p = —p. (5) 

Their solutions 

x(t) = xoexpit), p{t) = poexp{-t) (6) 

describe, for any E 0, the classical trajectory of a hyperbola E = Xopo- The system is 
unstable because it has a hyperbolic fixed point at x = 0,p = 0. 

The system is quantized by considering the dilation operator in the x space 

H = ^{xp + px) = -ih (^xd.^ + , (7) 

which is the simplest formally Hermitian operator corresponding to the classical Hamiltonian 
(4). 

The eigenvalue equation 

HiljEix) = E^e{x) (8) 

is satisfied by the eigenfunctions 

C 

where the complex constant C is arbitrary, since the solutions are not square-integrable. 
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The momentum eigenfunctions are found by the Fourier transform 

1 /■°° 

(f)Eip) = —7f iJE{x)exp{-ipx/h)dx. (10) 
V h J-oo 

To calculate the integral the function iPe{x) has to be analytically continued across the singu- 
larity at X = 0. Following [10], one might simply choose even eigenfunctions, which is a natural 
choice, since the classical hyperbolic Hamiltonian (4) is parity-symmetric. The result is (by 
using the reflection and duplication formulas for the gamma function): 

C p exp(-m) 

\p\ y z-K J~'x> \u\ ' ' 

Thus the quantum Hamiltonian (7) admits two independent stationary solutions [12] for 
any real energy E 

^e{±x) = C5(±x)e('^/'*-^/2)'°^l"l (12) 

together with the corresponding momentum eigenfunctions. In (12) 6{x) is the Heaviside step 
function. In fTD] a more general relation between coordinate and momentum eigenfunction is 
considered; this is interpreted as two phase-space currents, the x current flowing out from the 
origin and the p current flowing into the origin. 

The aim of our study is to obtain a discrete energy spectrum for the Hamiltonian (7) by 
imposing some quantization conditions. Quantization rules may arise in a natural way by 
the symmetries of the Hamiltonian (4). An obvious symmetry is the invariance under scale 
transformations: ^ 

x' = Xx, p' = —p (13) 

A 

with A positive. It follows from the classical equations of motion that the parameter of dilation 
A corresponds to evolution after time T = log A. We shall consider the scale invariance of 
the hyperbolic dynamical system as a part of a larger group of transformations and disscuss 
the Hamiltonian (4) from the point of view of conformal representations theory. Conformal 
symmetry in one-dimensional x space contains translations x' = x + a, dilations x' = Xx, and 
special conformal transformations. We note that the special conformal transformation 

/ X ~\~ cx , , 

X = ^ (14) 

1 + 2cx + c^x^ ^ ' 

can be written down as a superposition of two conformal inversions R (i.e. Weyl reflections), 

Rx = -- (15) 

X 

and a translation T^cX = x — c, 

x' = RT.cR- (16) 



3 



These transformations are generated by the differential operators 

P = —idx, D = —id — ixdx, 

K = -2idx-ix^d^, (17) 

which form a reahzation of the Lie algebra of the conformal group in one dimension, isomorphic 
to the Lie algebra of S0{2,1) ~ SL{2,R); d is a complex number, a conformal (or scale) 
dimension that labels the irreducible representations (IR). We shall briefly review the conformal 
group in one dimension S0{2, 1) ~ SL{2,R) and its principal series representation [T^, which 
will be used in the following. The group S0{2, 1) is the group of isometries of the pseudosphere 

-vo + yl + yl = -R\ 

the image of which is the unit disc, by a stereographic projection onto the plane. The conformal 
group in one dimension, known as the real Moebius group S0{2, 1) = SL{2, R)/Z2, acts as a 
group of automorphisms of the compactification of the real space TZ, which can be realized 
as a set of 'light' rays belonging to the cone Uq = yl + yi- This set of rays is in a one-to-one 
correspondence with the unit circle in two dimensions. The S0{2, 1) generators Xab = —Xba, 
a, b = 0,1,2 are expressed in terms of the conformal group generators by the relations: 

P = Xio + Xi2, K = Xu-Xio, D = X2o. (18) 

The image of the /^-inversion (15) on the cone is a reflection of the axes yi,y2, i-e. space 
reflection. In fact the stability subgroup of the point x = (or ?/i = 0,yo = y2 = ±1 in 
homogeneous coordinates) is isomorphic to the Weyl subgroup. The latter is a subgroup of 
SL{2, R) of order 2 generated {mod ± 1) by a non-trivial element w, with w'^ = —1. 

The unitary representations of G = SL{2, R) are induced by an unitary representation of its 
maximal compact subgroup K = 5*0(2) of rotations (generator X12). The representation theory 
is constructed in the space of functions with domain in the homogeneous space X = K \ G 
and values in a separable Hilbert space. The factor space X is isomorphic to the upper half- 
plane. It can be realized as the space of homogeneous functions of two variables of given parity 
and degree of homogeneity, determined by the conformal dimension —d. There is a one-to-one 
correspondence between the points of X and the points on the real line TZ^. An important 
class of unitary IR is the so-called principal series (and its dual) characterized by the value 
—d = —1/2 + ip (and —d = —1/2 — ip respectively), with p real. Let C°° be the space of 
infinitely different iable (complex-valued) functions f{x) of given parity, which satisfy certain 
asymptotic conditions for |a;| ^ 00 following from the smoothness of / in the neighbourhood 
of the i?-inversion element. The unitary (even) principal series IR acts in the Hilbert space 
completion Ti.^ of C°° with respect to the scalar product 

/oo _ 
Mx)f2{x)dx, (19) 
-00 

which is positive-definite and G-invariant. One can use any 5*0(2) invariant measure to con- 
struct a unitarily equivalent representation on the upper half-plane. 
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There exists a pair of intertwining maps relating the representation T^,x = "^P ^^'^ the dual 
one T*, X = —ip- It is defined by the convolution integral 



/■oo 



(20) 



J —oo 



and has the intertwining property 

The kernel of the integral operator (20) is given by 



y It: 




-l+2ip 



(21) 



The Fourier transform of G 



■oo 



GAx)e-'P''dx, 



(22) 



— oo 



restricts the choice of the normalization constant N^^ in such a way that, for equivalent repre- 
sentations, the product G^G^ is proportional to the unit operator 



corresponds to the Laplace - Beltrami operator L [14j, which is the second-order differential 
operator associated with the SL{2, i?)-invariant metric on the upper-half plane z = x + iy 
and has the form L = —y^A = —y'^{dl + dy). The solutions of the equation Lf = d{l — 
d)f are invariant automorphic forms with respect to the discrete subgroups F of SL{2,R). 
Automorphic forms are determined by the relations f{'yz) = f{z), for 7 G F, which serve 
as boundary conditions for the operator L. They define the Laplace - Beltrami operator as 
a self-adjoint operator acting in a Hilbert space of square integrable functions with respect 
to the invariant measure in which a unitary representation of the group G is realized. The 
procedure of extending L, which is a positive, semi-bounded and symmetric operator with a 
dense domain, to a self-adjoint operator is known [T3]. For unitary representations containing a 
stability vector with respect to the maximal compact subgroup, the spectral decomposition of L 
and the associated eigenfunctions - proper and generalized - are intimately connected with the 
discrete subgroup. According to the rigorous results on the L-theory if the discrete subgroup 
is such that the space G/T has a finite volume, the continuous spectrum entirely lies on the 
real half-axes 1/4 < d{l — d) < 00, while the point spectrum lies on the line < d{l — d) < 00. 
The point spectrum is the set of singular points of the kernel of the resolvant of L for which 
the invariant bilinear functional in the representation space is degenerate. The automorphic 
eigenfunctions satisfying Lf = d{l — d)f are generalized functions and form a complete set of 
functionals in the space of linear functionals on the unitary representation Hilbert space. 



gAp)gAp) = 1. 



(23) 



The Casimir operator of 5*0(2, 1) 



C2 
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The coordinate eigenfunctions iI>e{x) belong to the space (denoted hereafter Q*) of hnear 
functionals on the even unitary principal series IR with p — E/h. The position functions 
also belong to the Hilbert space extension Q* D H^. 

Returning to the eigenvalue equation (8), we may interpret it as a condition for a scale 
invariance of a conformal wave function of conformal dimension d = 1/2 — iE/h. If we postulate 
invariance of the wave function iIj{x) under the full conformal group, this will give a constant 
eigenfunction. It can be easily verified that in addition to scale invariance the wave function 
iIj{x) is invariant under conformal i?-inversion: 



Formula (24) defines the coordinate function as an automorphic form of weight d in the 
space of the IR principal series representation, invariant with respect to the discrete Weyl 
subgroup of SL(2, R), where its non-trivial element is chosen to be the conformal i?-inversion. 

(Automorphic is used here in the sense of invariant functions, associated with a stationary 
discrete subgroup of automorphisms of a manifold, which is the homogeneous space of a Lie 
group.) 

The scale- and i?-invariance properties of ip^x) suggest that the quantum wave function of 
a chaotic dynamical system (4) is a scale-invariant function, i.e. a homogeneous function of 
degree —d and an i?-inversion invariant automorphic function of weight d, where —d is the scale 
dimension that labels the unitary irreducible principal series representation and is hence 
equal to —1/2 + iE/h, with E real. The complex-conjugated wave function ■?/'(,x) corresponds 
to the dual representation with —d = —1/2 — iE/h. The above suggestion means precisely 
that the wave function of such a chaotic system is related to the kernel (21) of the intertwining 
operator for the IR principal series and its dual. As such it belongs to the space Jl* of linear 
functionals on the unitary principal series IR Hilbert space. This assumption assures that the 
eigenvalues of the formally Hermitian Hamiltonian (7) are real. 

We shall now show that the real values E are eigenvalues of a self-adjoint operator. For this 
purpose we consider the shifted Hamiltonian 1/2 — iH/h — —xdx, which satisfies 



The operator on the RHS of (25) is the Laplace - Beltrami operator L on the real x line. It is 
readily seen that the functions 



are eigenfunctions of this operator. We emphasize that the position eigenfunctions i(j{x) are even 
functions of x. As seen from the explicit expression (12) they are functions of the argument 
l^l only and can be considered as functions defined on the positive real line (0, oo). This 
assumption identifies ±x and allows the use of a parametrisation such that the representation 




(24) 




(25) 




(26) 
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space can also be realized as the Hilbert space of functions f{x), defined on the half line x > 
and square-integrable with respect to the measure x~'^dx. It is known [H] that the functions 

e{x, d) = + c(l/2 + ip)x^/^-'f, (27) 

where c{d)c{l — d) = 1, with the boundary condition at a with an additional parameter n, 

9{a) = aK9'{a), (28) 

form a complete set of eigenfunctionals of the self-adjoint operator L in the space of linear 
functionals on the unitary principal series IR Hilbert space of functions f{x) defined on the 
positive line a<x<cxD,a>l and square-integrable with respect to the measure x~'^dx. 

Let iP±e{^/x) be the eigenf unctions of the Laplace - Beltrami operator —x^dx corresponding 
to the eigenvalues 1/4 + E"^ /Tl^ of the continuous spectrum of L. It follows from (25) that 

xdx{l - x9,)a;i/2±is/^ = + |!^ ^i/2±iE/n_ (29) 

It is clear that x^^^^*^/^ are eigenvectors of the operators xdx and 1 — xdx^ corresponding 
to the eigenvalues 1/2 ± iE /fi and 1/2 =i= iE /fi respectively. Thus if 1/4 + E'^ /U belongs to 
the spectrum of L then 1/2 ±iE/h belongs to the spectrum of the dilation operator xdx and 
vice versa, so that roughly speaking the dilation operator is the square root of the Laplace - 
Beltrami operator in the Hilbert space of functions with domain (a, oo) and square integrable 
with respect to the measure x~'^dx. 

We turn now to the eigenvectors %l){x). As can be verified immediately they are eigenvectors 
of L but corresponding to complex eigenvalues 1/4 — (liiE'/^)^. However, there is an operator 
which, when acting on the position eigenfunctions '?/'(x), yields the real eigenvalues 1/ A + E"^ /fi^ . 
It is the i?-transformed Laplace - Beltrami operator denoted hereafter L^. It has the form 

Lr = -xdx{l + xdx) 

= -x^dl-2xdx. (30) 

Expressing Lr in terms of the Hamiltonian (7): 

Lr = {1/2 -iH/h){l/2 + iH/h) (31) 

we obtain that the position functions il){x) are eigenfunctions of Lr corresponding to the eigen- 
values l/A-\-E'^/1x^^ so that the operators l/2±iH/h are formally the square roots of a positive 
operator. The elliptic Laplace - Beltrami operator and the /^-transformed one act on functions 
whose arguments are related by the R transformation. The Laplace - Beltrami operator is 
densely defined and self-adjoint on the line [1, oo). The continuous /^-transformation x' = 1/x 
on the line x > 0, maps its eigenfunctions to the eigenfunctions of the /^-transformed operator 
Lr, which is hence densely defined and self-adjoint on the interval (0, 1]. The domains of the 
operators consist of elements satisfying the relation for automorphic /^-invariant functions and 
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have as common elements the functions / defined for x = 1, which is the fixed point of the trans- 
formation x' = 1/x. The two operators have to be equal on the common elements. Relations 
(24), which determine the eigenf unctions tplx) as invariant automorphic functions with respect 
to the conformal /^-inversion, serve as boundary conditions at the point a; = 1 (or at any point 
a > 1, if x' = a/x), for coincidence of the inversed operator Lji = {l — iH/h){l + iH/h) with the 
self-adjoint Laplace - Beltrami operator and will be used in phase space as quantization con- 
ditions that generate a discrete spectrum. Thus due to the ±x identification the i?-invariance 
automorphic relation defines the dilation Hamiltonian as a square-root of a self-adjoint oper- 
ator. For the eigenfunctions of the Laplace - Beltrami operator, eq. (24) is equivalent to the 
Neumann boundary conditions for the functions and for their first derivatives: 

(32) 

at the point |x| = 1. Cleary thanks to the additional paramater in (28) these conditions are 
trivially fulfilled. If not only the first of these conditions is automatically satisfied, while the 
second one has only the trivial solution, unless it is assumed that |x| = 1 is a point where 
orientation is reversed. This seems to be a natural assumption, since conformal /^-inversion is 
an orientation- reversing transformation, with fixed points |x| = 1, so that the points x may 
belong to manifolds of opposite orientation. Such are the curves of constant energy in phase 
space. The volume element in phase space is given by the area xp = h oi the Plank cell with 
sides Ix and Ip. Contours of constant energy in phase space are the upper and lower branches of 
the hyperbola E = px and of the conjugated one E = —px, corresponding to the eigenvalues of 
the of the square root operator ±H of the /^-transformed Laplace - Beltrami operator. Being 
aware of the time non-invariance of the classical orbits |10] , one should interpret the upper and 
lower branches of the conjugated hyperbola respectively as the combined space-time transform 
and the time transform of a given one. The /^-inversions in the x-space are easily enlarged 
to transformations in phase space, which are area preserving and leave invariant the operators 
±H. These transformations have the form: 

rr± , h , X^P 

T±: x' = -- P' = T^- (33) 

p h 

The upper it-indices denote that the determinant of the transformation is ±1. The — Ti trans- 
form of X is the conformal /^-inversion corresponding to the value a = 27r in (28) (in units 
h = 1). For xp = h the transformations contain the canonical one x' = —p,p' = x and re- 
fiections x' = ±x,p' = ±p. Together with the transformations —T^ and —T^ (note that the 
identity is —T2), they generate a finite group which is the dihedral group P4 of order 8 [15] . 
consisting of all orthogonal transformations in two dimensions that preserve the regular 4-sided 
polygon centered at the origin. It is actually generated by refiections (denoted hereafter Si and 
5*2 respectively) about the x axis and about the (diagonal) line / inclined at an angle 7r/4 to 
the positive x-axis. The rotations through multiples of 27r/4 form a cyclic subgroup of order 4. 
The invariant 4-gon is spanned by the vectors (±\/X, 0), (0, ±v^), (±v^, ±Vh). A fundamen- 
tal region in the square is the 2-simplex with vertices (0, 0), {\/h, 0), {\/h, \/h). The 2-simplices 
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are joint two by two to form four 2-complexes (the Plank cells in the distinct quadrants). The 
transformations fall into two classes depending on whether they preserve (reverse) sign of vol- 
ume element, i.e. preserve (reverse) orientation. We can use the isometries that map one side 
of the fundamental region onto a side of an adjacent region to glue distinct regions together 
along their boundaries, imposing boundary and quantization conditions and forming a closed 
path corresponding to a classical closed orbit. By identifying all boundary points mapped to 
one another by an element of the group (i.e. points in each orbit of the group) we can impose as 
a boundary condition the requirement that the coordinate and the momentum eigenfunctions 
are automorphic invariant forms of the discrete subgroup of reflections: 



for all s in the group. This is possible since for xp = ±/i, the x and p transformations each form 
a one dimensional factor map x' = T{x), p' = T{p) of the two dimensional map (33) of the 4- 
gon. The reflections in the x-factor map are the conformal i?-transformation and the a;-inversion 
and generate the Weyl subgroup of the dihedral group. The coordinate function ip{x) is thus 
defined to be an invariant automorphic function with respect to the discrete i?-transformation 
(which it is) and eq. (34) is equivalent to 



Eqs. (35) impose relations between values of and (p at boundary points mapped into each 
other by an element s, generating the finite Weyl subgroup of the dihedral group. Since under 
the discrete /^-transformation points transform to E > h with decreasing x and to E < —h 
with increasing x and taking the limit h ^ yields the semiclassical approximatiom, it will 
be thus imposed as a boundary condition for a smooth transition from E < —h to E > h. An 
orientation is then assigned to each 2-cell and a direction to each edge, the natural direction, 
due to the defining homeomorphism to the unit interval [0, 1], being from initial point to final 
point 1. To each distinct element gi,i = 1,2, ..8 (i.e. root vector) of the group a vertex i is 
assigned and a directed edge connects two vertices iff gi = SkQj for a generator Sk of the 
group. A path extends from a given g at an edge of the fundamental simplex to Si-^...Sif,g and 
is closed if it can be presented in the form S'j^...S'jj, = 1. This is a consequence of the relation 
(SiSj)"^^'''^^ = 1, with m{i,i) = 1, for the Coxeter element [16] of the dihedral group, which 
is just the product of two genetating reflections, hence is a rotation through 27r/4 of order 4. 
Denoting the rotation through 7r/2 by r = 5*1 5*2, we obtain that a closed path S'l 5*2 5*1 5*2 is always 
of the form rr, which is the word representation of a projective plane. Thus the relation for the 
Coxeter element forces identification of antipodal boundary points and introduces topology of 
a projective plane [I7| in phase space. At both 1^. and Ip boundary sides of the 4-gon ±x and 
±p are identified, the latter according to the equivalence relation {—\/h,p) ^ {Vh, —p), which 
endows the phase space with Moebius topology. It is readily seen that the Weyl reflection 5*2 
and S2 form a cyclic subgroup of order 2, isomorphic to the group Z2. With the identification of 
boundary points a closed path begins from a given (initial) point gi = {x, p) and ends at a (final) 
antipodal one gj = (— —p). The closed paths 7^ fall into two homotopy classes depending 
on whether the final point is reached via a product of even or odd number of reflections 82- In 



ip{x) = ip{sx) 



(j){j>) = (t){sp) 



(34) 





(35) 
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the first case, the end point is connected to the initial one by a transformation with det — 1 
and orientation is preserved along the path, while in the second case, the determinant of the 
transformation is —1 and orientation is reversed along the path. Hence the two homotopy 
classes form a representation of the fundamental group of the projective plane HiilZV'^), which 
is the non-trivial homology group Z2. The boundary conditions for the automorphic functions 
have to be compatible with the topology of the projective plane, and in particular with the 
Moebius topology equivalence relation. Due to the twist the original 2-complexes (the Plank 
cells in the original quadrants) are mixed and eq. (35) becomes 

'^i^)lv=h = (-^) \{-x)i-p)=h = ^V' (^) U-p)=h (36) 

The boundary conditions in the form of eq. (36) are compatible with the identification of ±a;, 
and of ±p by a twist, and serve at the same time as quantization conditions generating discrete 
spectrum. Indeed we have: 



C 



^l-2iE/h ^-l/2+iE/?i 



Using further that x — h/ ± p and the expression (11) for the momentum eigenf unctions, we 
obtain 

4^(x) - ±^.V2..^^^.i.A£(lZi^l^ .38) 



With the help of the functional equation for the Riemann zeta function 

,_.,,r((i-^ 

sv J r{s/2) ' 



(39) 



for s = 1/2 + iE /h, we get from (38) the result 

a;V2^(l/2 _ iE)'il}{x) = ±p^/^C(l/2 + iE/h)(l){p), x{±p) = h. (40) 

The latter formulas, with ±x and ±p identified, are consistent only if either ^jj{x) = (and 
hence (f){p) = 0), or 

C(l/2 + iE/h) = C(l/2 - iE/h) = 0, (41) 

which provides a regularization procedure for semiclassical approximation and generates a dis- 
crete real spectrum of the dilation Hamiltonian operator (7). Geometrically the conditions 
refiect the fact that the fundamental group of the projective plane is non-trivial, and can be 
mapped homomorphically to the automorphisms group of order two of the covering space. The 
defining homomorphism 

a : 7^l(7^P^x) ^ ±1 ~ Z2 

assigns to each path 7^ a number (7(75) = ±1, according to whether orientation is preserved or 
reversed by transport around it in such a way, that (j(7s) 7^ 1, for 7^ 7^ 1. This is consistent 
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with the group structure of the transformations (33), which form the double point dihedral 
group with the factorization property V4 = T>2 x Z2 (it can be presented equivalently as a 
semidirect product of the two cyclic groups Z4 ■ Z2)[T8]. Hence the irreducible representations 
have a Z2-grading too. They are twice as many, and given a representation i(j{x) of V2 each 
element of the pair cr(7s)'?/'(x) could be interpreted as associated with an orientation preserving 
or reversing homotopy class. For x belonging to a closed trajectory the quantization conditions 
(40) are an algebraic expression of the homomorphism a, this unique homomorphism being the 
geometric representation of the finite Weyl (orthogonal) reflection group. 

Employing the notation as in [10] 

(xlV^)^V^(x), im^m, H = H/h, (42) 

we rewrite expression (40) in the form: 

{x\x'/\{l/2 - iH)\ijE) = ±W^Cill2 + iH)\ijE). 

xp = h, {—Vh,p) ^ {Vh, —p)- (43) 

These relations can be geometrically interpreted as quantization conditions generating the 
Riemann zeros as a discrete spectrum of the energy operator of the conformally invariant 
Riemannian metric due to dihedral symmetry in phase space of the chaotic dynamical system. 
The Riemann zeros are the eigenvalues of the hyperbolic chaotic Hamiltonian for which the 
latter is the square root of the (inversed) Laplace - Beltrami operator. 

The relation with the — sign on the RHS of (43) is exactly the quantization condition 
suggested by Berry and Keating in (TU] . 

Remark: It is known that the surface symbol is not unique. Hence we can interpret the 
closed path either as a cylinder 6*1 6*2 5*2, which is obtained from the 4-gon by identification 
of the Zj^-sides only, or as a Moebius band S1S2S1S2 with the /p-sides identified by a twist. In the 
first case the boundary condition is the trivial one (the + sign in (43)) and is not a quantization 
condition. In the case of a Moebius band the boundary condition is the non-trivial one (the — 
sign in (43)) and serves as a quantizatin condition which generates Riemann zeros. 

To summarize we have implemented conformal symmetry to model a chaotic dynamical 
system. The quantum Hamiltonian has a classical limit and a discrete real spectrum whose 
points are the Riemann zeros on the critical line Res = 1/2. The coordinate eigenfunctions 
are automorphic functions, invariant with respect to the discrete subgroup of Weyl reflections, 
for points on hyperbolic trajectories of constant energy in phase space. The smooth semiclas- 
sical behaviour around the saddle point in phase space is only possible for eigenvalues of the 
Hamiltonian, which are the critical zeroes of the Riemann zeta function. 
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